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A.  PURPOSE  OP  PROQRAM 

To  compute  the  loss  coefficient  of  a  reservoir  outlet  works  for  a  given  gate 
opening  from  drawdown  data  given  the  stage  versus  reservoir  volume  relation 
and  knowing  the  average  inflow  rate  (may  be  zero)  during  the  observation  period 
of  drawdown. 

Reference:  Moore,  E.  T. ,  "Reservoir  Storage  Evacuation  Time  by  Integration," 
Civil  Engineering,  April  1971- 

6£Jr'  ' 


C.  METHOD*  'B..  ' 

The  program  is  written  in  G6 35  time-share  series,  FORTRAN  IV,  and  is  part  of 
a  Conversationally  Oriented  Real-Time  Program-Generating  System  (CORPS).  The 
program  consists  of  a  main  program.  This  main  program  handles  all  I/O  re¬ 
quirements,  computations,  and  writing. 


D.  EQUIPMENT  DETAIL* 

The  program  was  developed  and  is  operational  on  the  WES  G635,  Vicksburg,  MS. 
It  is  also  operational  on  HIS  66/80,  Macon,  GA,  and  Boeing  CDC,  Seattle,  WA. 


Input  requirements  are:  observation  period  of  drawdown  in  hours,  initial  and 
final  gage  heights  of  water  surface  in  ft  above  outlet  intake  invert,  outlet 
conduit  cross-sectional  area  in  ft2,  coefficients  of  power  curve  to  data  of 
gage  height  (ft)  versus  reservoir  storage  (acre-ft),  and  the  average  inflow 
into  the  reservoir  in  acre-ft/hr.  Output  includes  the  given  input  data  and 
the  loss  coefficient  of  the  outlet  conduit. 


P.  ADDITIONAL  REMARK* 

Complete  documentation  of  this  program  is  available  from  the  Engineer  Computer 
Programs  Library,  Technical  Information  Center.  WES. 
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B.  PROGRAM  SPECIFICATIONS: 

Language :  ANSI  FORTRAN  (FORTRAN  IV) 
Solution  Requirements:  The  run  command 


RUN  WESLIB/CORPS/H2251,R 


and  the  inputs  defined  in  E. 

Method  of  Analysis:  Numerical  integration  of  the  continuity  equa¬ 
tion  combined  with  the  stage  versus  outflow  and  stage  versus  reser¬ 
voir  storage  relationships. 

Core  Requirements  G635:  11  K  words 

External  Storage:  None 

Restrictions :  Period  of  observation  of  drawdown  should  be  less  than 
24  hours  to  meet  computational  limitation  of  the  GE  637  computer 
system. 

General  Equation:  Continuity  equation 


where  Q  is  the  outflow,  B  is  a  constant  inflow  and  —  is  the 
change  in  storage. 


REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appertdix  B 

PART  I:  ENGINEERING  DESCRIPTION 

1.  PROGRAM  NUMBER:  722-F3-R0-2D1 

2.  TITLE :  H2251  -  Computation  of  Loss  Coefficient  for  an  Outlet  Works 

from  Prototype  Drawdown  Data. 

3.  REVISION  LOG:  N/A 

U.  PURPOSE  OF  PROGRAM:  To  compute  the  loss  coefficient  of  a  reservoir 
outlet  works  for  a  given  gate  opening  from  drawdown  data  given  the  stage 
versus  reservoir  volume  relation  and  knowing  the  average  inflow  rate 
(may  be  zero)  during  the  observation  period  of  drawdown. 

Reference:  Moore,  E.  T. ,  "Reservoir  Storage  Evacuation  Time  by  Integra¬ 
tion,"  Civil  Engineering,  April  1971. 

5.  STEP  SOLUTION: 

a.  Input .  The  input  variables  that  must  be  defined  prior  to  exe¬ 
cuting  the  program  are: 

(1)  The  observation  period  of  drawdown  in  hours,  i.e.,  the  time 
required  for  the  reservoir  water  surface  to  fall  from  an 
upper  level  to  some  lower  level.  The  time  should  be  limited 
to  less  than  2h  hours  to  avoid  computational  errors  de¬ 
scribed  in  paragraph  5c. 

(2)  The  initial  and  final  gage  height  of  the  reservoir  water 
surface  in  ft  above  a  datum  such  as  the  invert  of  the  out¬ 
let  works  intake.  Any  datum  is  satisfactory  provided  it 
is  also  used  in  the  gage  height  versus  reservoir  volume 
relationship  in  (U)  below. 
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(3)  The  cross-sectional  area  of  the  outlet  conduit  in  ft  . 


(4)  The  coefficients  of  a  power  curve  fit  to  the  plot  of  gage 
height  (ft)  versus  reservoir  volume  (acre -ft).  A  straight 
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line  is  fitted  to  the  log-log  plot  and  the  equation  has 
the  form: 


Y  =  AX8 


where  Y  is  the  reservoir  volume  (acre-ft),  X  is  the 
gage  height  (ft)  (the  same  datum  use  in  (2)  above  is  re¬ 
quired),  and  A  and  B  are  the  coefficients  required  by 
the  program. 

(5)  The  average  inflow  into  the  reservoir  during  the  observa¬ 
tion  period  in  acre-ft  per  hour.  Zero  inflow  is  the  most 
desirable  inflow  condition  because  it  simplifies  the 
solution. 

b.  Theoretical  Formulation.  Formulation  of  the  problem  is  derived 
from  the  continuity  principle:  i.e.,  outflow  minus  inflow  equal  to  a 
change  in  storage  or: 


(l) 


where  Q  is  the  outflow,  B  is  a  constant  inflow  and  —  is  the 

at 

change  in  storage  or  volume  in  the  reservoir.  The  outflow  (Q)  is  de¬ 
fined  by: 

Q  =  A^^(cfs)  =  0.0826A.J^F  (acre-ft /hr)  (2) 

2 

where  A  is  the  cross-sectional  area  of  the  outlet  conduit  (ft  ),  h 

2 

is  the  total  head  (ft),  g  is  the  acceleration  of  gravity  (32.2  ft/sec  ) 
and  K  is  the  outlet  works  loss  coefficient.  Flow  rates  in  this  formu¬ 
lation  are  expressed  in  acre-ft/hr  because  of  the  common  practice  of  ex¬ 
pressing  reservoir  storage  in  units  of  acre-ft. 
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The  reservoir  volume  or  storage  is  expressed  by  a  power  function  of 
gage  height,  i.e.. 


V  =  Ch 


where  V  is  the  reservoir  volume  (acre-ft)  and  h  is  the  gage  height 
or  head  (must  correspond  to  h  in  equation  2).  The  coefficients  C 
and  D  are  determined  from  a  straight  line  fit  to  the  data  of  volume 
versus  gage  height  plotted  on  a  log-log  plot.  If  a  simple  straight  line 
does  not  fit  the  data  over  the  range  of  gage  heights,  straight  line  seg¬ 
ments  of  best  fit  should  be  used.  Thus,  C  and  D  of  the  segment  con¬ 
taining  the  values  of  observed  gage  height  would  be  selected. 

The  rate  of  change  of  V  with  respect  to  time  can  be  written  as 


dV  _  __  D-l  dh 
dt  -  CDh  d¥ 


Substituting  equations  2-U  into  equation  1  and  solving  the  resulting 
integral  equation  using  a  power  series  expansion  (see  Ref  page  l), 
an  expression  for  the  time  (hrs)  required  to  draw  the  reservoir  water 
surface  elevation  from  one  level  to  a  lower  level  with  constant  inflow 


is  obtained: 


Id) 


f2D-i  2D-i 

-i2 


2D  -  i 


where,  E  =  0.0826a^^&  ,  h^  and  h^,  are  the  initial  and  final  water 
surface  gage  heights,  respectively,  and  N  is  the  number  of  terms  in 


1  , 


the  series  expansion  (magnitude  of  N  depends  upon  desired  accuracy). 
If  the  constant  inflow  (B)  equals  zero,  then  the  terms  for  i  =  2,  3, 
...»  N  equal  zero,  thus  greatly  simplifying  the  solution,  i.e.. 


t 


(6) 


Substituting  the  definition  of  E  into  equation  >,  the  functional  rela¬ 
tionship  in  terms  of  the  outlet  works  loss  coefficient  (K)  is  expressed 
by: 


where : 


F(K)  = 


2CD 

0. 0826A v^ 


i  /  2D-i 
Zi  K2  (hT2 


t 


(7) 


Zi  = 


_ B1  1 _ 

(0. 0826A»^2g)1_1  (2D-i) 


(8) 


Note:  is  set  to  1/2D-1  if  inflow  (B)  equals  zero. 

c.  Numerical  Solution.  The  solution  of  equation  7  is  obtained  by  a 
first  order,  Newton  fixed  point  iterative  scheme 


F(K  ) 

Ki+1  =  Ki  "  FOO 


(9) 


where  is  the  new  approximation  of  the  loss  coefficient  based  on 

a  previous  approximation  (K^)  and  the  value  of  the  function  (equation  7) 
and  its  first  derivative  (F'(K^))  with  respect  to  K  .  Differentiating 


equation  7  with  respect  to  K 


H2251 


F'(K) 


CD 


0.0826A/2g 


2D-i 


Li=l 


(10) 


equation  9  as  two  terminating  criteria.  The  first  is  defined  by 


K..,  -  K. 
l+l  i 


<  e 


(11) 


where  e  is  a  user  specified  error  tolerance.  Its  value  in  the  program 

is  set  at  0.01.  The  second  cirteria  is  that  if  convergence  specified  by 

* 

equation  11  is  not  obtained  in  200  iterations,  the  program  stops.  A 
third  terminating  criteria  is  used  in  equation  7  to  minimize  the  number 
of  terms  (N)  in  the  series  since  the  GE637  is  limited  to  a  maximum  ex¬ 
ponent  of  38.  Computation  is  terminated  when  the  difference  between 
F(K)i_i  and  FCK)^  does  not  exceed  one  percent,  i.e.. 


<  0.01  (12) 


However,  if  equation  12  is  not  satisfied  after  35  terms  (i  =  35) »  the 
program  stops  and  a  message  is  printed  stating  that  "reasonable  accuracy 
cannot  be  obtained."  If  the  period  of  observation  is  restricted  to  less 
than  2U  hours,  the  criteria  in  equation  12  should  be  obtainable  with 
i  <  35. 


F(K).  -  F(K)._1 
F(K).  +  F(K)._x 
2 


d. 


Steps  in  Program  Execution. 

(1)  Input  variables  defined. 

(2)  The  initial  value  of  the 
to  1.0. 
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loss  coefficient  (K)  is  set  equal 


e>  w  ,  c*  e  a  as  *» 

>■  H-  •<:  r:  a  >h  o 

-j>  1  t.  a’  m  1-1  •  o 

r*  I  ,  r*  fi  Q  M  * 

n  b  -in 

1-  1  ••  !  ►♦on  w 

• '  r.  y*  o 
'  •  :  '  1  o 


y-' 


o 

o 
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(3)  If  the  inflow  (B)  is  zero;  is  set  equal  to  unity  and 

equations  7  and  10  are  solved  with  N  =  1  ;  go  to  step  5. 

(h)  If  the  inflow  (B)  is  not  zero,  equations  7  and  10  are  solved 
with  i  =  1  ,  2,  ....N  until  the  difference  between  F(K)^ 
and  f(K)jy_^  is  less  than  or  equal  to  1  percent.  If  N 
becomes  greater  than  35,  the  program  stops  and  an  error 
message  is  printed,  "reasonable  accuracy  cannot  be  obtained 
with  given  data." 

(5)  Newton's  method  is  applied  to  compute  another  approximation 
for  K  using  equation  9  and  the  values  of  F(K)  and  F'(K) 
from  step  3  or  step  h.  If  the  difference  between  the  new 
approximation  of  K  and  the  previous  value  of  K  is  less 
than  0.01  (equation  11),  the  solution  is  complete  and  the 
results  are  printed.  If  the  criteria  of  equation  11  is  not 
satisfied,  the  program  returns  to  step  3  or  step  ^  (depend¬ 
ing  upon  the  value  of  inflow)  with  the  new  approximation  of 
K  and  repeats  until  equation  11  is  satisfied  or  200  itera¬ 
tions  are  reached.  If  the  solution  requires  more  than  200 
iterations  the  program  stops  and  a  solution  is  not  possible 
for  the  given  input  data.  Once  equation  11  is  satisfied 
the  results  are  printed. 
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PART  II:  COMPUTER  FUNCTIONAL  DESCRIPTION 

1.  REVISION  LOG:  N/A 

2.  FUNCTIONAL  FLOW  CHART: 


ms 
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3.  EQUIPMENT  AND  OPERATING  SYSTEM:  The  program  was  developed  on  a  G635 
time-share  system  in  which  input/output  equipment  consisted  of  a  Model  3? 
remote  teletype.  It  is  now  operational  on  the  WES  G635»  Vicksburg,  MS; 

HI.'  i  '/'s'/,  Macon,  GA;  and  Boeing  CDC ,  Seattle,  WA. 

1.  MM  Ml  IHKMKNTS :  Tne  required  inputs  are  entered  via  the  user's 
time-share  terminal  device  in  free-field  format. 

'  •  SECONDARY  SiUKAGK  INPUT:  None. 

v.  MMM'T  DATA  DESCHI1TI0N:  The  following  names  are  used  for  the  input 
variable;-  in  program  H2251. 

2 

.•IKEA  -  cross-sectional  area  of  outlet  conduit,  ft1- 
COEKKA  -  y  intercept  coefficient  of  power  fit  curve 
CuEFEB  -  slope  of  power  fit  curve 
FLOWIN  -  average  inflow,  af/lir 
HEAD1  -  initial  gage  height,  ft 
HEADS  -  final  gage  height,  ft 
TIME  -  period  of  observation,  hr 

7.  OUTPUT  DATA  DESCRIPTION:  The  follow  name  is  used  for  the  output 
variable  in  program  H2251. 

K  -  loss  coefficient 

8.  PROGRAM  MESSAGES:  Return  from  all  program  messages  is  RERUN  or  STOP . 

a.  Reasonable  accuracy  cannot  be  obtained  for  data.  Suggest  de¬ 
creasing  time  of  observation  in  data  collection. 

b.  Newton  method  exceeds  200  iterations  -  stopped. 

,  t 

i 


o 
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VARIABLE  DEFINITIONS: 

A  -  temporary  variable 

A A  -  temporary  variable 

AAI  -  temporary  variable 

2 

AREA  -  cross-sectional  area  of  outlet  conduit,  ft 

COEFFA  -  y  intercept  coefficient  of  the  power  fit  curve  in  the  gage 
height  (ft)  versus  reservoir  volume  (af)  relation 

COEFFB  -  slope  of  power  fit  curve  to  the  gage  height  (ft)  versus 
reservoir  storage  (af)  relation 

DIFF  -  percent  difference  (fraction)  between  the  value  of  F(K)  , 
equation  7,  and  F(K).. 

DK  -  difference  between  the  new  approximation  to  the  loss  coef¬ 
ficient  (K)  and  the  previous  value  of  K 

E  -  temporary  variable 

EXP1  -  exponent,  i  -  1 

EXP2  -  exponent,  2D  -  i/2 

EXP3  -  exponent,  i/2 

FK  -  functional  relation  in  terms  of  K  ,  equation  7,  with  the 

number  of  terms  (N)  in  the  summation  equal  to  i 

FK1  -  same  as  FK  except  the  number  of  terms  (N)  in  the  summation 

equal  i  -  1 

FKP  -  the  value  of  the  first  derivative  of  F(K)  with  respect  to 

K  ,  equation  10 

FLOWIN  -  constant  inflow,  af/hr 

HEAD1  -  initial  gage  height  of  reservoir  water  surface,  ft 

HEAD2  -  final  gage  height  of  reservoir  water  surface,  ft 

HFILE  -  5  character  name  of  program  (H2050);  passed  to  WESLIB 
routine  HACCT  for  bookkeeping 
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JKL  -  directs  return  from  WESLIB  routine  RERUN  to  desired  input 
read 

K  -  outlet  works  coefficient 

KSTAR  -  new  approximation  to  loss  coefficient  K  ,  equation  9 

LQX  -  equal  1,  print  instructions  from  WESLIB  routine  RERUN; 
equal  3»  no  print 

LQZ  -  equal  1,  execute  all  input  cues  and  reads;  equal  2,  call 
WESLIB  routine  RERUN  and  execute  only  desired  reads 

SUM1  -  value  of  the  terms  summed  in  equation  7  for  the  initial 
gage  height  and  N  equals  I  -  1 

SUM11  -  same  as  SUM1  except  N  equals  I 

SUM2  -  value  of  the  terms  summed  in  equation  7  for  the  final  gage 
height  and  N  equals  I  -  1 

SUM22  -  same  as  SUM2  except  N  equals  I 

SUM3  -  value  of  the  terms  summed  in  equation  10  and  N  equals  I 
TIME  -  observation  period,  hr 

ZI  -  value  of  equation  8  for  each  increment  of  i 

ZZZZ7.  -  2  character;  equal  to  RE,  rerun;  equal  ST,  stop 

iO.  EXAMPLE  CASE: 

Input:  Observation  period  (TIME)  =  U.l67  hr 
Initial  gage  height  (HEAD1 )  =  71  ft 
Final  gage  height  (HEAD2)  =  70  ft  ^ 

Cross-sectional  area  (AREA)  =  16  ft 
Coefficients  of  power  curve  (COEFFA)  =  .0^35 

(COEFFB)  =  2.57^ 

Average  inflow  (FLOWIN)  =  0  af/hr 


OUTPUT : 
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H225 1  -  CALCULATION  OF  0"TLET  VOPX?  LOSS  COEFFICIENT 
FPOM  PPOtotvoe  MEAS”FEMENT  OF  DFAWDOVN 

AA-EV-EP  "EPIOD  OF  OBSERVATION  <HPS>. 

-4.167 

AP-EN'-EP  INITIAL  GAGE  HEIGHT-FEET  ABOVE  OUTLET  INVEPT. 
-71 

AC-EMTEP  FINAL  GAGE  HEIGHT-FEET  ASO"E  OUTLET  INVEPT. 
-70 


AD-ENTEP  CPOSS- SECTIONAL  APEA  OF  OUTLET  CONDUIT  IN  FT**2. 

-16 

AE-EWTEP  COEFFICIENTS  OF  PO'-’EP  C"pv£  FIT  TO  GAGE  HEIGHT!  FT)  V«  PESEPVOIP 
tOL»lNE<  ACpE-FT)  .  FOPN  OF  EQUATION:  Y-A*X**&>"  WHEPE  v  EOUALS  GAGE  HEIGHT. 
X  EOUALS  VOLUME  AND  A  AND  B  APE  THE  COEFFICIENTS.  ENTEP  A  AND 
THEN  P  SEPEPATED  VI TH  A  COMMA. 

-.0435,2.574 

AF-ENTEP  THE.AVEPAGE  INFLOV  PE®  HOUP  IN  ACPE-FT/HP  FOP  THE 
PEPIOD  OF  OBSEpuflTION. 

-0 


INPUT 

PEPIOD  OF  OBSEFUATI ON  -  4.17  HPS. 

INITIAL  GAGE  HEIGHTCFT.  above  OUTLET  INVEPT)  “  71.00 

FINAL  gage  HEIGHT! ft.  ABOVE  Outlet  INVEPT)  -  70.00 

CPOSS- SECTIONAL  APEA  OF  OUTLET  CONDUIT  -  16.00  FT--2 

COEFFICIENTS  OF  POVE®  CPVE-GAGE  HEIGHT  VS  PESEPVOIP  'VOLUME. 
PO  PM  OF  EQUATION  *  Y-A-X—B  VHEPE: 

A  “  0.04350 

B  -  2.57400 

INFLOV  -  0.  ACPE  FT/HP 


OUTPUT  H2251J 

LOSS  COEFF.  -  16.7128 


ENTE®  PEP"N  OP  STOP 
■STOP 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  C 

PART  III:  FILE  DOCUMENTATION 

1.  REVISION  LOG:  N/A 

2.  TITLE :  H2251  -  Calculation  of  Outlet  Works  Loss  Coefficient  from 

Prototype  Measurement  of  Drawdown. 

3.  PROGRAM  SOURCE  LISTINGS:  See  pages  lU-17. 

1*.  NUMERICAL  AND  LOGICAL  ANALYSIS: 

Newtons  first  order  fixed  point  iteration  scheme  is  used. 

5.  SUBROUTINES  NOT  DOCUMENTED  IN  ABSTRACT:  None 

6.  MISCELLANEOUS :  The  program  is  part  of  the  CORPS  computer  system. 
CORPS  is  an  acronym  standing  for  Conversationally  Oriented  Real-Time 
Program-Generating  System.  The  program  is  now  operational  on  the  WES 
G635,  Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC ,  Seattle,  WA. 
The  source  listing  on  page  lU  contains  the  first  line  run  command  and 
brief  for  H2251.  This  first  line  run  command  runs  the  binary  H2251B  of 
the  source  listing  on  pages  15-17  (FORTRAN  source  of  H2251)  and 
attaches  to  WESLIB  routines  HACCT  and  RERUN. 
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03/02/79 

H225I 


000  -.»E'LI  P/'l22S  1  P,  r>}  "E'LI  B/ CEP"*I,  o;  VE'L  I  P/H  ACC',  P 

0800  69  rn  P  P*>0GcAM  COMOit-nre  -HE  LOP'  COF"I  Cl  E*J'  OF  A  pEeEr,,'OIP  0 ”T- 
08  1  0  65LF'  '.’0 'K '  F0B  A  Gl*'E'J  GA'E  OPEMIVG  F«*)M  rBAVVE0V'J  DATA/  G1  "EN  THE 
08  20  6  9<?-AGE  .t-  eECjmf)  j  o  'fOL'ClE  PEL  ATI  0*1  AHD  KMOWPIG  'H  E  A'^AGE  IV- 
0830  69FL0”  "A'E  <*1A*'  PE  TEPO)  D"PIMG  ‘'HE  0EFEP,'A'I0M  PEPIOD  OF  DPAW- 
0390  05pO<«j. 

08  5  0  6  2  imojitc  DEO”IPED  APE:  DPAUMDOVH  0P'EP"A'I0*I  PEPIOD  HP,  INITIAL 
08  6  0  6  2a*ID  FIHAL  GAGE  HEIGH"  OF  '-'A'EP  '"'FACE  F'  AP0"E  0'"LE'  IM'AKE 
03  7  0  61 1  *I*»E",  0,,'LE'  C0*ID'»I'  G°0  Pf-  PEC'10‘1  A'EA  'OF*,  COEFFI  Cl  EV'P  OF 
03  3  0  6  2pO"E”  C”'"E  TO  DATA  OF  GAGE  HEIGH'  FT  VP  PEFEP”OI«>  "OPAGE  AF, 
0890  93 AND  -HE  A"EPAGE  IHFLOV  IM'O  'HE  PEPEP"0IB  AF/HP. 

09  0  0  6  9  0"T°"'  I?ICL"DE?  THE  GI"EM  DA'A  AHD  THE  LO"  COEFFI  Cl  EH'  OF  THE 

09  10  I  SO'f'LE'  C0*JD"IT. 

09  9  9  *  06fI*1I  P*I 
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0000 1  *  |t!"M  *  =  J  /C9rr,'/H225  lBCIOGO) 

10000  Lf7=  liL™*  UKKK=6J  PEAL  K,KFTAP 
10010  CHA0ACTE0  HFILE*5;  H*I L E» 5HH225 1 
10020  ppt'ir  7 

10030  7  P0fMA*’(//"H225l  -  CALC”LA-I3'J  OF  0'”LET  VOPKF  LO"  COEFFICIENT 
100404"/  1  Ik,  ••FPOM  PPOTOTVPE  MEA?*»PEMEMT  OF  DPAVEOVN"/) 

10050  15G00  CALL  :iACC"(HFILE) 

10060  GO  TO<  5,3  5) ,  LOT 
10070  5  ppIM*  10 

10030  10  FOPMA'TC/,,AA-EMTEP  PEPIOD  OF  OEFEPVATIOM  CHPF>*"> 

10090  PEAD,  TI ME 
10100  GO  "OC 1 5,85) ,LOT 
10110  15  ^PI«I~  20 

10120  20  *-OeMA~C/,'AS-E7!,,E?  INITIAL  GAGE  HEI GHT-  FEET  ABOVE  0 VTLET  Pl'lEPT. 
101  304") 

10140  FEAD,  HEAD  1 
10150  GO  TO  <  25,  8  5)  ,  LOT 
10160  25  30 

10170  30  FOPMAT<  /"AC-EMTEP  F1MAL  GAGE  HEIGHT-FEET  ABOVE  0"tlEt  INVEPT.") 
10130  °EAD«  HEAD2 
10190  GO  T3<05,85),LCT 
10200  35  o»I'lT  40 

10210  40  FO PMAt< /"AD-ENTEp  CPO 5 F- SECTI OHAL  APEA  OF  OUTLET  C0ND"1T  IN  FT* 
102204*2.-) 

10230  "FAD. APEA 
10240  GO  TO  <  45,  3  5)  ,  LOT 
10250  45  ppiMT  50 

10260  50  FOFHATC /"AE-ENTE*  COEFFI CIENT'  OF  "OVEP  CVP”E  FI  T  ~0  GAGE  HEIGH 
102704T(FT)  "t  pEFEpvo  I  P"/"”OL’F1E<  ACPE-  F*)  .  FO*M  OF  EC-AnON:  ','*A*X**PJ 
102304VHEPE  v  EC'AL '  GAGE  H El  GH T,  "/"K  EOVALF  "OL’PIE  AMD  A  AMD  B  AFE  THE 
102904COEFFICI  E’|tp.  ENTEP  A  AMD  ••/••'HEM  E  PEPEPATED  WITH  A  COMMA.") 

10300  PEAD, COEFFA, COEFFE 
10310  GO  TO<  55, 85), LOT 
10320  55  pr>i»iT  60 

10330  60  FO°MAt( /"AF-57JTEP  THE  A"E?AGE  INFLOW  PEP  HO"P  IM  ACPE-FT/HP  FOP 
103404  THE'V"”EPIOD  OF  OBFEPVATI OM.  ”) 

10350  PEAD, FLOVIM 

10360  GO  TO  C  9  5,  8  5)  ,  L^T 

10370  35  CALL  PE?"MCKHK,Lor» JKL) 

10380  GO  t3(  5,  I  5,  25,  35,  45,  55,  95),  JKL 

10390  95  nopjT  100,  tt  *JE,  HEAPl »  H  EAP2,  A  PEA,  COEFFA,  COEFFP,  FLO  VIM 
10400  100  FO"MAT(  /  "I  MPVT"//  "pepI  0  P  .  0  F  OFFEP"A.tIOM  F7.  2,  "H®F.  "/ 

1041 04"I MI ti AL  GAGE  HEIGHt(FT.  AP0"F  0"”LE-  IMVEPT)  *  ", F7. 2, /, "FIMAL  GA 
104204GE  HEIGHTCFT.  APO”E  0"TLEt  IM”Ept)  »  ",  f7.  2, /,  "C®0  **■- rEC  TI  OVAL  APE 
104304A  0*  0"”LEt  COND"!*  *  ", F7.2, "  Ft** 2", /,  "CO EFFI  C I  EH  TP  OF  POVEP  C"” 
104404 UF- GAGE  HEIGHT  VF  PEPEp-OIP  T»OL"ME.  "/,  ”F0  °M  OF  ECA'ICM:  v*A*',»*E 
104504"HEpps  "/,  5y,  "A  ■  ",  Fl  0.  5, /,  5x,  ”E  *  ",  Fl  0.  5,  /,  "IMFLO  V  »  ",*7.2,"  AC 
104604nE  FT/HP"/) 

10470  E»COEFFB-0.S 

10430  AA"0. 08  264*APEA*FO”T< 64.  4) 
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H225  1  c  COMT. 


10490  flAI *COEFFA*C0EFF3/AA 

13500  K=<  <TME«E)/C  AAI*(HEADU*E-HEAD2**E)  ) >**2. 

10510  I«-(<n.O'>lM.LE.0.0)  GO  TO  130 
10520  JJ*0 

10530  8  A=AA/ «ApT(K) 

10543  c*Tl  1*0.  0»  «pm2*0. 01  <tTi3*0.  Ci  I  ■  1 
10SS0  9  Exr>l  =  rLOA'r(  I- 1 ) 

10560  ETn2=  t  2.  *  CO  SFFB-  FLO  ATf  I  > ) /2« 

10570  Ev"3»FL0AT<n/2. 

10530  »I  =  (n.O”!*l**BXPl  )/<  AA**£XM*2.*r<P2) 

10590  et/Ml  l  =  F'*Ml 

10600  jmm 1  =  cmm l  +  lI*K**tE'CPl/2.)*HEADl**E'tP2 
10610  5»"l22=<;"M2 

10620  ?tiM2»e:itM2+r  i*k**<  ETPl/2. )  *H  EAD2«*  EM  P2 

10630  F'»M3*«"M3+lI*F1.0ATt  1  >*X**<  FLOATt  I)/2.-l . 0)*<HEACl**  EMP2-HEAD2«*EXp 
10640*2) 

10650  FKP*AAI*P"H3 

10660  FKlc  t  2.*C0EFFA*C0EFFS*t  ?<fMl  1-  FPM22)  J/A-’TME 
10670  FK=  t  2. *COEFFA*  COF.FFE*  t  SUM  1  -  S!’M2)  )  /  A-T1ME 
10680  IFtI.EO.l)  GO  TO  116 
10690  mrF*CFK-FKl)/C<FK+FKl)/2.> 

10700  ircABs<mrr).LE.0.0i)  go  to  120 

10710  116  1=1+1 

10720  I  Ft  I  • GE. 35)  p»IMT  121 

10730  lFd.GE.3S5  GO  TO  135 

10740  121  FORMAT! /"PEA 50MASLE  ACCtTACY .  CAWOT  BE  OBTAINED  FOP  DATA.”/” 
107504,«''GGE«:t  DECFEA«IMG  TIME  OF  OBFEP’'ATIOM  IN  DATA  COULEC^IOM.  "/> 

10760  GO  TO  9 
10770  120  CO*l*-IM*'E 
10730  KFTAP=X-FK/FKP 
10790  DK”ABc<K-K?TAP) 

10300  JJ“JJ+l;  KaK FTAP 
10310  1F< JJ. GE. 200)  GO  TO  125 
10320  IFtDK.GE. 0.001)  GO  TO  8 
10330  GO  TO  130 
10340  125  PPITIT  126 

10350  126  FO®MATC/"MEVrO'J  METHOD  EXCEEDED  200  I  TEPATIOMS-FH)PPED"/> 

10360  130  PFIMT  1 53, K 

10370  150  FO°MA-C//"OPTPt'T  H2251  !'V”LO?p  COEFF.  »  ",Fl0.4) 

10380  135  LCT=2 

10890  CHAPACTF.p  t’.TTE*2 

10900  16000  PPI'JT  16002 

10910  16002  FOPMATC//"E',JTEP  PEPITJ  OP  STOP") 

10920  PEAB  16001,  TTTTT 
10930  16001  F0PMAT<a2) 

10940  I  Ft  TTTTE  •  EG.  2'IPE)  GO  T0  15000 
10950  I  Ft  TT7.5T  •  EG.  2H  FT)  GO  TO  20000 
10960  PFlWr,  "EPPO®  ***  PE*"'PE” 

10970  GO  TO  16000 
10980  20000  CTOP 
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